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Abstract 

We present a next-to-leading order computation in QCD of one-jet and two-jet cross 
sections in polarized hadronic collisions. Our results are obtained in the framework of a 
general formalism that deals with soft and collinear singularities using the subtraction 
method. We construct a Monte Carlo program that generates events at the partonic 
level. We use this code to give phenomenological predictions for pp collisions at y/S = 
500 GeV, relevant for the spin physics program at RHIC. The possibility of using jet 
data to constrain the poorly known polarized parton densities is examined. 


^Work supported by the Swiss National Foundation. 




1 Introduction 


In the last few years, measurements |l[] of the spin asymmetries in deep-inelastic scat¬ 
tering (DIS) of longitudinally polarized lepton beams off polarized hadron {N = p, n, d) 
targets have provided much new information on the spin structure of the nucleon. Theoreti¬ 
cal leading order (LO) [2-6] and next-to-leading order (NLO) [2-7] analyses of the data sets 
demonstrate, however, that these are not sufficient to accurately extract the spin-dependent 
quark (Ag = — q^) and gluon (Ag = — g^) densities of the nucleon. This is true 

in particular for Ag{x,Q^), since this quantity contributes to DIS in LO only via the Q^- 
dependence of , which could not be accurately studied experimentally so far. As a result 
of this, it turns out [2-6] that the a:-shape of Ag seems to be hardly constrained at all by the 
DIS data, even though a tendency towards a fairly large positive total gluon polarization, 
/q Ag{x, = 4 GeV'^)dx > 1, was found Biia- The measurement of Ag thus remains 
one of the most interesting challenges for future high-energy experiments with polarized nu¬ 
cleons. In selecting suitable processes for a determination of Ag, it is crucial to pick those 
that, unlike DIS, have a direct gluonic contribution already at the lowest order. Here, one 
thinks in the hrst place of high-pj. reactions in nucleon-nucleon collisions, which have been 
tremendously important in the unpolarized case to constrain the unpolarized gluon density. 

At the moment, the most eagerly awaited experimental tool for the ‘spin physics’ commu¬ 
nity is the RHIC collider at BNL, at which hrst runs in a proton-proton mode with polarized 
beams are expected to be accomplished in about two years from now. The centre-of-mass 
energy for these pp collisions will be ranging between 100 and 500 GeV, with luminosities 
(rising with energy) between 240 and 800 pb“\ respectively. One expects about 70% po¬ 
larization for each beam. Such conditions look extremely favourable for studying the spin 
asymmetries for all kinds of high-p^ pp processes that are sensitive to the gluon density, such 
as jet, prompt-photon, or heavy-havour production. Considering the higher end of RHIC 
energies, jets could be the key to Ag: at y/S = 500 GeV, clearly structured jets will be ex¬ 
tremely copiously produced, and jet-observables will show a strong sensitivity to Ag thanks 
to the dominance [^ of the gg and qg initiated subprocesses in some kinematical ranges. 

In order to make reliable quantitative predictions for a high-energy process, it is crucial 
to determine the NLO QCD corrections to the Born approximation. Quite in general, the 
key issue here is to check the perturbative stability of the process considered, i.e. to examine 
the extent to which the NLO corrections affect the cross sections and (in spin physics) the 
spin asymmetries relevant for experimental measurements. Only if the corrections are under 
control can a process that shows good sensitivity to, say, Ag at the lowest order be regarded 
as a genuine probe of the polarized gluon distribution and be reliably used to extract it from 
future data. 

NLO QCD corrections are expected to be particularly important for the case of jet- 
production, since it is only at NLO that the QCD structure of the jet starts to play a role in 
the theoretical description, providing for the hrst time the possibility to realistically match 
the experimental conditions imposed to dehne a jet. 

The calculation of the NLO QCD corrections to jet production by polarized hadrons 
is the purpose of this paper. Such a calculation needs the one-loop 2 —2 and tree-level 
2 —3 polarized (i.e. not summed over external helicities) amplitudes as input. Fortunately, 
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these amplitudes are already known |^, |^. Furthermore, several independent methods to 
calculate any infrared-safe quantity in any kind of hard unpolarized collision are at present 
available in the literature BB P!3| . The formalism of ref. m has been used in ref. [Q to 
construct a Monte Carlo code that can calculate any three-parton infrared-safe observable 
in hadron-hadron unpolarized collisions. In the present paper, we will extend the method 


of refs. [^, Q and adapt the Monte Carlo code to the case of polarized hadron-hadron 
collisions. As a result, we will present a customized code, with which it will be possible to 
calculate any infrared-safe quantity corresponding to either single- or di-jet production to 
NLO accuracy. 

The outline of the paper is as follows: In section ^ we describe the calculation of next- 
to-leading order corrections to jet cross sections in polarized hadronic collisions. Since there 
is no conceptual difference to the unpolarized case, we restrain from repeating all details 
and only give an overview. A crucial input for the next-to-leading order calculation is 
the polarized gluon distribution function. Section ^ summarizes the current situation and 
discusses the main assumptions that are used to constrain the various hts. We begin our 
phenomenological study in section with the investigation of the perturbative stability and 
reliability of our next-to-leading order results for single-inclusive as well as double-differential 
observables. In section § we turn to the more phenomenological issue of the dependence of 
the cross sections on the various parametrizations of the parton densities. In particular, we 
will also investigate some spin asymmetries. Finally, in section]^ we present our conclusions. 
Some technical details concerning the difference between the calculations of the cross section 
in polarized and unpolarized hadronic collisions are given in an appendix. 


2 Formalism 


We start by writing a generic differential jet cross section in polarized hadronic collisions 


using the factorization theorem |15 




aia2 A1A2 


X daaia 2 {xiKi, \j; X 2 K 2 , \ 2 ;S), 


( 2 . 1 ) 


where Hi and H 2 are the incoming hadrons, with momenta Ki and K 2 and helicities Ai 
and A 2 respectively; is the non-calculable but universal distribution function for the 

parton a, with helicity Aj in the hadron Hi with helicity Aj and dd'a^a 2 ci-re the (subtracted) 
short-distance partonic cross sections^]. The quantity S is the measurement function; it 
dehnes the jet momenta in terms of the parton momenta. Its specihc form therefore depends 
upon the jet-clustering algorithm adopted. Since the results for polarized scattering are 

^In the following, we will often omit some of the entries of the quantity ddaia 2 J th® meaning should be 
clear from the context. 
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usually presented in terms of the quantity 




( 2 . 2 ) 


it is convenient to rewrite eq. (^) in the following form: 


K2-, S) = '^ dxidx2Af^^^\xi)Af^^^\x2)dAaaia2{xiKi;x2K2;S), 


aia2 


where dAaaia 2 is defined analogously to dAa^^^^'^^ in eq. (E3). and 

A f(Hi) _ AHi+) _ p{Hi+) _ AHi—) _ AHi—) 

ai j Cli-\- j CLi— JcLi— j CLi-\- 


(2.3) 

(2.4) 


The available parametrizations for polarized parton densities are presented in terms of the 
quantities Afa- 

Our aim in the present paper is to evaluate the polarized cross section in eq. (p.3|) at next- 
to-leading order in perturbative QCD, for one- and two-jet production. As is well known, 
there are two main problems in such a computation. First of all, the relevant polarized am¬ 
plitudes have to be known. For our case, we need the 2 —>■ 2 one-loop and 2 —>■ 3 tree-level 
amplitudes. These amplitudes have been computed P, |T^ using the technique of colour 
ordering and the helicity method. Secondly, the cancellation of the infrared poles, which 
appear in the intermediate steps of the calculation, has to be performed analytically. Due to 
the universal (i.e. process-independent) structure of these poles, general methods |TT, II3.0 
exist, which allow the computation of any infrared-safe cross section in any type of hard scat¬ 


tering. In the present paper we will use the approach of ref. , based upon the subtraction 
method. The main idea of ref. |]T3 is to exploit the properties of the measurement function 


in order to disentangle the infrared-singular regions that appear in the real contribution. 
Indeed, as follows from the universal properties of the measurement function (which are 
responsible for the good dehnition of any infrared-safe cross section), this quantity can be 
written as a sum of terms, each of which is non-vanishing only when one given parton is 
soft or collinear to another parton. Therefore, the partonic cross section can be expressed 
as a sum of terms whose singular structure is trivial, and the subtraction procedure can 
be straightforwardly implemented. It is important to notice that, although the approach of 

H 


ref. 


was originally introduced in the case of unpolarized collisions, it does not need any 


principle modification to be applied to the case of polarized collisions, since the major role 
in the treatment of the singularities is played by the measurement function. Further details 
on this topic can be found in the appendix. 

The procedure of ref. 


12] results in the subtracted partonic cross sections which appear 


in eqs. G) or 
quantity 


Standard Monte Carlo methods can therefore be used to compute the 


{H)= / dAa^^^^^\S) H, 


(2.5) 
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where H is any function of the jet momenta, which are dehned by S. If H is a product of 
9 functions, implementing experimental cuts and selecting a bin of a given histogram, then 
{H) is the QCD prediction for the cross section in that bin. 

The main drawback of eq. (|2.5|) is that the jet definition is used, through the measurement 
function, to disentangle the infrared singularities. This prevents us from getting predictions 
for different jet definitions when performing a single computer run, as is customary in a par- 
ton shower Monte Carlo; also, one-jet and two-jet observables must be treated separately. 
This problem was addressed in ref. [H; the key observation is that, in order to disentangle 


the infrared singularities, the measurement function can be substituted by a suitable sum 
of products of 9 functions, which we call the V function. By construction, each term of 
the sum is non-vanishing only in one given infrared-singular region which contributes to the 
cross section at next-to-leading order. In this sense, the V function and the S function are 
completely equivalent, and the former can be interpreted as a fake measurement function 
(however, while the real measurement function contains 5 functions that define the jet mo¬ 
menta in terms of the parton momenta, the V function only contains 9 functions). Therefore, 
following ref. 


14| , we can write the analogue of eq. 


as 


IffS] = j 


( 2 . 6 ) 


By construction, we get {H) = [i^5]. A computer program based upon eq. (|2.5|) outputs 
jet momenta, which are eventually used to £11 a histogram as specified by H. On the other 
hand, a computer program based upon eq. (|2.6| ) outputs parton momenta; these quantities 
are eventually used to compute jet momenta (as specified by S), which will again fill the 
histogram given by H. Therefore, in eq. ( ^.bj) no jet definition is involved in the generation 
of the hard event and in the computation of the weight. It follows that, during the same 
computer run, it is possible to use several different jet definitions and to compute one-jet, 
two-jet and non-jet (like transverse thrust) observables. This kind of computer code is called 
a parton generator. 

two parton generators were presented, one for photon-hadron collisions and 


In ref. 


one for hadron-hadron collisions. The latter has been suitably modified to deal with polar¬ 
ized hadron-hadron collisions, and used to produce the phenomenological results presented 
in this paper. The structure of the code of ref. [1^ remains unchanged, since only the 


unpolarized partonic cross sections and splitting functions had to be substituted with the 
polarized ones (see the appendix for more details). 


3 Polarized Parton Distribution Functions 

As stated in the introduction, there is hardly any experimental information on the spin- 
dependent gluon density Ag at present. In contrast, the quark densities are far better 
constrained by the existing data from inclusive polarized DIS. This is particularly true for the 
polarized valence quark densities, which come out rather similar in all theoretical analyses 
performed so far. The spin-dependent sea-quark distributions seem less well constrained; 
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Figure 1: The polarized gluon (left) and valence quark densities (right), as given 
by the six NLO parametrizations that will be used in this paper, at the scale = 
100 GeV^. The patterns for the quark densities are the same as those used for the 
gluon. 


however, they are of minor importance for our jet studies^. 

In our phenomenological analysis we will try to cover as much as possible of the wide 
range of polarized parton densities allowed by the present DIS data: this is especially relevant 
for the gluon density, as we will discuss below. We will use the following six polarized parton 
density sets, which were obtained within various theoretical analyses of polarized DIS: 

• Two sets of ref. |0, the ‘standard’ set, corresponding to the best £t to the data obtained 
in ref. [|| (from now on, referred to as GRSV std), and a set obtained when saturating 
the positivity constraint lA^fj < g for the gluon density at the low input scale 0(0.5 
GeV) (GRSV maxg). 

• The three sets of ref. [§] (DSSl, DSS2, DSS3), obtained by constraining the hrst 
moment of the polarized gluon densities in three different ways. The hrst moment of 
the gluon of DSSl is the largest and almost one order of magnitude larger than for 
DSS3, which has the smallest gluon. 

• Set G of ref. P] (GS-G), which provides a gluon distribution with a qualitatively dif¬ 
ferent x-shape, becoming negative at large x for low Q^. Sets A and B of ref. ^ are 
similar to some of the GRSV and DSS ones, and we do not use them so as to avoid a 
proliferation of curves. 

^Strictly speaking, inclusive DIS data for proton and neutron targets can only give information on the two 
non-singlet and the quark singlet combinations, rather than on all quark and antiquark densities individually. 
The inclusion of data ^ for the spin asymmetry in semi-inclusive DIS (SIDIS), which in principle would 
allow a complete flavour discrimination in the quark sector, does not provide much help in practice, since the 
SIDIS data have not yet reached the precision of the inclusive ones. The flavour decomposition of the nucleon 
parton densities - which is essential for making predictions for observables other than DIS - therefore partly 
depends on theoretical assumptions made when analysing the DIS data [2-5]. 
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Figure 2: Ratio of NLO to LO polarized gluon densities. We use the same pattern 
as in fig. to distinguish the various lines. Also shown is the ratio for an unpolarized 
(GRV) set. 


All these distributions are available at both LO and NLO, the latter corresponding to the 
MS scheme used also in our calculation of the NLO partonic cross section. 

In this work, we will analyse the phenomenology of polarized jet production for the 
particular case of RHIC at a centre-of-mass energy of y/S = 500 GeV and for jet transverse 
momenta in the region of 15 GeV < < 100 GeV. This implies that the polarized parton 

distributions will mainly be probed in the x-range 0.06 < x < 0.4 and at typical scales 

of the order of 100 GeV^. Figure shows the NLO polarized valence quark and gluon 
densities of the six different sets we are going to use. It becomes obvious that there is indeed 
a wide range of possible gluon distributions compatible with present polarized DIS data. 
As mentioned earlier, in the valence sector, most of the distributions are very similar, with 
slight exceptions in the cases of the Uy distribution of GS-G and the dy density of GRSV 
std. The origin of these differences can be easily traced back to the fact that the analyses of 
refs. [ii are based on a somewhat smaller data sample than the others, since some data 
sets became available only after [^, ^ had been published. The spin-dependent sea-quark 
densities (not shown in £g. differ more strongly among the various sets, but have only 
a very small impact on the jet cross sections in this kinematical region: they only amount 
to less than 5% of the contribution from valence quarks. In conclusion, since the variations 
in the quark sector are much smaller than the ones for gluons, we can expect that any 
differences between predictions for the polarized jet cross sections (or asymmetries) that 
are found when using different polarized parton density sets, are to be attributed to the 
sensitivity of the observable to Ag. 

The size of radiative QGD corrections to a given unpolarized hadronic process is often 
displayed in terms of a ‘iF-factor’ which represents the ratio of the NLO over LO results. 
In the calculation of the numerator of K one obviously has to use NLO-evolved parton 
densities. As far as the denominator is concerned, a natural definition requires the use 
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of LO-evolved parton densities. However, by using NLO-evolved parton densities and LO 
partonic cross sections, one still obtains a hadronic cross section accurate to LO, and therefore 
the denominator of the iL-factor can also be computed with NLO-evolved parton densities. 
If one chooses a ‘natural’ subtraction scheme, such as MS, these two dehnitions of the K- 
factor are expected to give similar results (we stress that the two dehnitions might give 
rather different results in the framework of an arbitrary subtraction scheme: there is no 
reason to worry about that, since the iL-factor is not a physical quantity, and therefore it 
is not supposed to be scheme-independent). However, in the case of polarized scattering, 
additional problems arise. Indeed, suppose one attempts to ht Ag from the DIS data. 
Since the data hardly constrain the gluon density, very different results for Ag can emerge 
if the ht is performed at LO or at NLO. This is conhrmed in hg. where we show the 
‘gluonic’ iL-factors = Ag^^^/Ag^^ as functions of x for our various set^. It can 
be seen that indeed for most sets K^g is not close to unity. We observe that things are 
much better in the unpolarized case, where there are far more data to constrain the gluon: 
here Kg = ^nlo^^lo ^ Also note that in the case of the ‘GRSV maxg’ set, where one 
assumes that = ^lo.nlo input scale (see above) in both LO and NLO, 

a Kj\g 1 is found also at higher scales. This underlines our point that artihcially large or 
small iL-factors for, say, polarized jet production could result merely from the fact that the 
gluon is at present so ill-constrained. We will discuss this point further in section |^. 


4 Perturbative Stability 


As we mentioned in the previous section, in the case of polarized collisions the study of 
the iL-factor, which in general does not give any information on the perturbative stability 
of the results in hadronic scattering, faces additional problems, due to the large difference 
between and To illustrate this issue, we have computed the LO hadronic 

cross section entering the iL-factor in two different ways, by convoluting the partonic LO 
cross sections with either LO-evolved parton densities or NLO-evolved parton densities. For 
the sole purpose of distinguishing the two dehnitions, we will call the former the tree-level 
cross section and the latter the Born cross section. As an example we have chosen the 
spectrum, using the Ellis-Soper (ES) cluster jet algorithm as proposed in ref. with the 
resolution parameter D = 1. In hg. ^ we show the ratio of the NLO cross section over the 
tree-level cross section. This ratio can be rather large for some parton density sets. However, 
a comparison with hg. ^ shows that these large corrections come mainly from the change of 
Ag^^ to Indeed, if we plot the ratio of the NLO cross section over the Born cross 

section, as done in hg. ^b, we see that the corrections are moderate. For the same reason as 
in hg. D we do not show the curve for the GS-C set. On the other hand, we also show the 
same ratio for the unpolarized case, using the GRV parton densities . 

A reliable error estimate on our NLO results requires some knowledge on the size of the 
uncalculated higher-order terms. The only fully reliable way to get this information is to 
perform computations of even higher order. Unfortunately, such a calculation is currently 
out of reach. Thus, the best we can do is to study the dependence of the full NLO results on 


^We do not show for GS-C, since in this case may be zero. 
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Pt (GeV) Pt (GeV) 

Figure 3: Ratio of the next-to-leading order cross section over (a) the tree-level 
cross section (i.e. leading order pdf) and (b) the Born cross section (i.e. next-to- 
leading order pdf) for various parton densities. We use the same pattern as in fig. 
to distinguish the various parton densities. 


the renormalization and factorization scales. Throughout we will set the two scales equal, 
i.e. fiR = fip = fx. Although physical observables are obviously independent of /i, theoretical 
predictions do have such a dependence. It arises from the truncation of the perturbative 
expansion at a hxed order in the strong coupling constant Og. A large dependence on p, 
therefore, implies a large theoretical uncertainty. 

The scale dependence of the tree-level and Born cross section is very similar. For the rest 
of this paper we will always use the Born cross section (i.e. the leading order hard partonic 
cross section integrated over the next-to-leading order parton densities) as the leading order 
result. As we will see, the scale dependence is substantially reduced once the next-to-leading 
order corrections are included. 

We will always consider proton-proton scattering with y/S = 500 GeV. For the strong 
coupling constant Og we use the standard two-loop form with Aqcd set to the value used in 
the parton distribution function under consideration. Our default choice for the scale is 



i 


(4.1) 


where the sum is over all hnal-state partons and kij. is the transverse momentum of parton i. 
We will study the scale dependence of several observables by comparing the results obtained 
with with those obtained with /i = /io/2 and /i = 2/io, in both the polarized 

and unpolarized cases. For this purpose, we £x the parton distribution functions: our 
default choices are MRST |]T3 and GRSV std P for the unpolarized and the polarized case 


respectively. We verihed that our conclusions are unchanged when using other density sets. 









Pt (GeV) Pt (GeV) 


Figure 4: Scale dependence of the next-to-leading order and Born pr-distributions 
for the Ellis-Soper algorithm with D = 1. (a) Polarized pp scattering and (b) unpo¬ 
larized pp scattering at y/S = 500 GeV. The range of the pseudo-rapidity is restricted 
to \r]\ < 1. 


To start with, we consider the inclusive p^-distributions, where Pt is the transverse 
momentum of the jet. As in £g. we use the cluster jet algorithm with the resolution 
parameter D = 1 and require \ri\ < 1. We checked that, basically, the dehnition of a jet 
through a cone algorithm only amounts to a change in the normalization. Concerning the 
scale dependence, the different jet dehnitions give very similar results. In £g. ^ we show the 
next-to-leading and leading order distributions in the polarized case for the three different 
scales: po, 2/io and Note that the Born results have been rescaled in order to disentangle 
the two sets of curves. Figure^ is the corresponding plot for the unpolarized case. Clearly, 
the dependence on the scale is substantially reduced when going to next-to-leading order. 
The situation in the polarized case is indeed very similar to the unpolarized one. 

We also investigated the distribution in pseudo-rapidity and verihed that the scale de¬ 
pendence is reduced in a similar manner as for the pr-distribution. Again, the dependence 
on the scale is very similar to the unpolarized case. 

A more stringent test on the perturbative stability of hxed-order QCD calculations can 
be made by considering quantities that are more exclusive than single-inclusive ones. At 
next-to-leading order, this basically means double-differential cross sections. Although the 
definition of these quantities is to some extent arbitrary, we will use the most commonly 
adopted prescription: we select all the events with two or more jets, we apply suitable cuts 
to the two leading jets (i.e. those with the largest transverse momentum) of each event, and 
we hnally compute correlations between these two leading jets. To be specihc, we require 

Pit > PTt^ Ihll < 1, P2t > P2T > Ih2| < 1, (4.2) 

where Pi^ and pi denote the transverse momenta and the pseudo-rapidities of the two leading 
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Figure 5: As in fig. for the xi-distribution (see the text for the definition). 


jets. We considered the two cases 


p'ff) = pfft = 10 GeV, 


(4.3) 


and 


^cut 


Plfi = 10 GeV, 


p™* = 15 GeV, 


(4.4) 


and studied various correlations, such as the angular distance in the plane perpendicular to 
the beam axis the invariant mass Mjj and the transverse momentum pfi! of the pair, 

the difference in pseudo-rapidity of the two jets /S.pjj and so on. In the case when the cuts 


given in eq. (^) are applied, we hnd that the next-to-leading order results are more stable 
(with respect to scale variations) than the corresponding Born ones, where this comparison 
makes sense (there are regions of the phase space that are not accessible at the Born level. A 
typical example, with the cuts of eq. (E3). is the threshold in the invariant mass distribution. 
In these regions, the NLO results display a scale dependence that is larger than anywhere 
else. This is what we expect, since in these regions our results are effectively leading order 
results). As an example, we present in £g. ^ the distribution in the variable Xi, dehned as 
follows 

Pit -t-P2T 

xi =- -= -. (4.5) 


The variable xi roughly corresponds to the momentum fraction of one of the two partons 
of the scattering process, and it is exactly so at the Born level. From the hgure, we see 
that the NLO result has an enhanced scale dependence in the region of small Xi; this is 
due to the fact that the small Xi values correspond to almost back-to-back jets, that is to 
a conhguration which is sensitive to soft gluon emission. However, in the whole remaining 
range of xi the scale dependence of the NLO result is smaller than that of the Born result. 

The case when the cuts given in eq. (|4.3|) are applied is more problematic in the frame¬ 
work of a hxed-order QGD calculation. This situation has been discussed in great detail in 
ref. 1^; here, we just remind the reader that, although the cuts given in eq. (|4.3|) define 
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an infrared-safe cross section (and can therefore be implemented withont any problem by 
the experiments), there are regions in the phase space (basically, all those corresponding to 
an exactly back-to-back conhgnration of the two leading jets) where any hxed-order QCD 
calcnlation breaks down, and an all-order resnmmation wonld be reqnired. We explicitly 
verified that, in the regions A(pjj — vr, pl^ ~ 0 and close to the threshold in Mjj, the scale 
dependence of the NLO resnlt is mnch larger than that of the LO resnlt, thus signalling a 
failure in the perturbative expansion. However, we stress that, in all the remaining regions 
of the phase space, the cuts of eq. (|4.3|) result in a cross section as well behaved as that 
obtained by applying the cuts of eq. i3)- 

Finally we would like to mention that we have computed all the aforementioned distri¬ 
butions by using both the ES and the cone jet-hnding algorithms. We observed only minor 
differences between the two resulting cross sections. 


5 Dependence of Observables upon Parton Densities 


Up to this point, we have shown that the polarized jet cross sections at next-to-leading 
order display a remarkable stability under scale changes in the upper end of the energy 
range probed at RHIC. It is therefore sensible to use our code to investigate in some detail 
a few phenomenological issues relevant to hadronic physics at RHIC. 

We thus turn to the problem of studying the dependence of the theoretical predictions 


upon the perturbatively non-calculable quantities that enter eq. (|2.1|) , namely the parton 


densities and the value of Aq^o- We will not perform an analysis of the separate dependence 
of our results upon these two quantities, since the polarized parton density parametrizations 
are only available with a single value of Aqcd- In the following, Aqcd will therefore always 
be set equal to the value associated to the various parton density sets used. In this section, 
our predictions are obtained by defining the jets with the ES algorithm with D = 1, and the 
scales have been set to the default value p. = /iq. 

We start by considering the asymmetry 


dAa/dpT 

da/dpr 


(5.1) 


Here Aa and a are the one-jet inclusive cross sections for polarized and unpolarized scattering 
respectively, and p^ is the transverse momentum of the observed jet. A cut |p| < 1 has also 
been applied. In Eg. ||, Apj, is shown as a function of Pt- The results for Aa have been 
obtained by choosing the six different parametrizations of the polarized parton densities 
previously mentioned. The unpolarized cross section a has been evaluated using the MRST 
set. Figure ^ clearly shows that the choice of the polarized parton densities induces an 
uncertainty on the theoretical results of more than two orders of magnitude. This enormous 
spread is basically due to the fact that the polarized gluon density is very poorly constrained 
by present DIS data, and at this energy the jet cross section is dominated by gg- and 
ggf-initiated parton processes. Therefore, there is a chance that the measurement of the 
polarized jet cross section at RHIC will be useful in order to rule out some of the polarized 
sets that are at present consistent with the data. However, since we see from Eg. 0 that the 
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Pt (GeV) 


Figure 6: Asymmetry versus transverse momentum for various polarized parton 
densities. The minimum measurable value for the asymmetry is also shown. 


asymmetry is always rather small, regardless of the specihc densities used, the statistics will 
have to be large. To estimate the minimum value of the asymmetry observable at RHIC, we 
use the well-known formula 




1 1 


(5.2) 


where C is the integrated luminosity, P is the polarization of the beam, and the factor e < 1 
accounts for experimental efficiencies; a is the unpolarized cross section integrated over a 
small range in transverse momentum bin). The quantity dehned in eq. (|5.2|) is plotted 
(boxes) in £g. |], for e = 1, P = 0.7, C = 100 pb“^ and apy-bin size of 2 GeV. From the hgure 
we see that the asymmetry defined in eq. (|5.1|) is measurable if the polarized parton densities 
are as described by most of the sets considered here. On the other hand, if the densities are 
as suggested by the GS-C set, the measurement will be possible only by combining a large 
statistics with a good overall efficiency (clearly, it is also possible to enlarge the bin size, in 
this way decreasing the minimum observable asymmetry, at the price of losing resolution). 
As a general feature, we may notice that the slope of the minimum observable asymmetry is 
steeper than the slope of the theoretically predicted asymmetries, increasing with increasing 
transverse momentum. This means that, although the value of the asymmetry is larger at 
high Pt than at small Pt, the measurement in this region will be more problematic. We also 
point out that, apart from the case of GS-G, the shape of the asymmetries obtained with the 
different parton density sets is rather similar, which just reflects the fact that the various 
gluon densities are rather similar in shape. 

We can still argue that the values for the asymmetry displayed in fig. are artihcially 
small, because of the different values of used when computing Aa and a. In the latter 
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Figure 7: Asymmetry versus pseudo-rapidity for various polarized parton densities. 
The minimum measurable value for the asymmetry is also shown. 


case, the MRST set has = 220 MeV, while in the former case the typical value is around 
= 130 MeV. For this reason, in the inset of £g. we plot the ratio of the unpolarized 
cross sections obtained with the MRST and the GRV sets (the GRV set has = 131 MeV). 
We see that if we had plotted the asymmetries using GRV instead of MRST, the final result 
would differ by 15% at the most, with respect to what is shown in fig. Our previous 
conclusions are therefore unaffected by the choice of the unpolarized densities. Finally, we 
verified that, if one defines the jets using a cone algorithm instead of the ES algorithm, 
the results are practically unchanged: as one might have anticipated, the details of the jet 
definition ‘cancel’ in the ratio that defines the asymmetry. 

The definition given in eq. (|5.1|) can be generalized to any single-inclusive or double¬ 
differential observable. In fig. we plot the asymmetry for single-inclusive jet production, 
as a function of the pseudo-rapidity of the jet, with a cut Pt > 15 GeV. As should be 
clear from the small-pr region of fig. also in this case all the predictions obtained with 
different density sets lie above the minimum observable asymmetry calculated with e = 1 
and C = 100 pb“^ (the bin size is 0.2). As in the case of transverse momentum, all the 
shapes are rather similar, except for the one relevant to GS-G, which has a local minimum 
a.t 7] = 0. The study of the asymmetry in the region around p = 0 could therefore be used to 
infer information on the shape of the polarized gluon density. We verified again that fig. ^ 
is unchanged if we use a jet-hnding algorithm based upon a cone prescription. 

We also computed single-inclusive asymmetries at the Born level. The results differ 
from those presented in fig. ^ and fig. for a factor up to 20%. The shape is also different. 
Therefore, NLO corrections give non-trivial information on the structure of the asymmetries. 
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Figure 8: Asymmetry versus transverse momentum. The effect of a cut on the 
azimuthal distance between the two leading jets of the event is shown. 


The possibility remains that, in some regions of the phase space, the asymmetry is larger 
than that shown in hgs. ^ and ^ and therefore that the measurement of the polarized cross 
section turns out to be easier. In the following we give a simple example. By inspection of 
our results for double-differential observables, we noticed that the correlation is steeper 
in the case of unpolarized collisions than in the case of polarized ones. In other words, the 
fraction of events with Acpjj ~ tt is higher in unpolarized collisions with respect to polarized 
collisions. This means that, by rejecting all the events with Acfjj ~ tt, the asymmetry will 
be enhanced. We therefore considered all the events with at least two jets, and we evaluated 
the azimuthal distance Acfjj between the two leading jets. We eventually rejected all the 
events with Acfjj > 3.05. For all the jets in the events passing this cut, we then computed 
the asymmetry as given in eq. (|5.1|) . The results are shown in £g. |], for two different parton 
density sets, superimposed on the results already shown in fig. obtained with the same 
parton densities. As we expected, the value of the asymmetry is indeed increased, by a 
factor of about 1.5 to 2. However, if we compute the minimum observable asymmetry as 
indicated before, we see that this quantity also increased (since the unpolarized cross section 
has decreased by effect of the Acfjj cut). It follows that the measurement of the asymmetry 
in the case of a Acfjj cut will be as difficult as in the case when no cut is applied (and, in the 
high-py region, even more difficult). On the other hand, such a possibility is still interesting 
in the sense that it can be used as a valuable test of the predictions of perturbative QCD, 
since the observable shown in fig. is more exclusive than the one presented in fig. ^ 

We finally consider some double differential observables, namely Acpjj, Mjj, Arjjj and the 
Xi-distribution, defined in the previous section; the results are presented in terms of Act in 
fig. p. For the plots shown in that figure we choose the cuts as given in eqs. (14.2|) and (4.4). 
Each observable is computed with the six different parton density sets used previously in 
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Figure 9: Double differential observables in polarized collisions for various parton 
densities. The order of the cross sections is the same as in fig |^. 


the case of asymmetries. We see that, for Acfjj, xi and rjjj, the shape of the distributions is 
rather insensitive to the choice of the densities. A somewhat more pronounced dependence 
is displayed in the case of the invariant mass Mjj. In fact, the typical value for the Bjorken 
X probed in two-jet events is of the order of We explicitly verihed that the pattern 

shown in £g. ^ is not modified if we dehne the jets with a cone algorithm. The same is also 
true if we change the scales from their default values, except for the Acfjj ~ vr region, where 
the scale dependence is very large and the shape of the distribution is sizeably modihed. 
This is due to the well-known fact that this region is particularly sensitive to the emission 
of quasi-soft gluons, and a resummation to all orders would be needed to get a reliable 
theoretical prediction. As discussed in section § we do not have a similar problem in the 
small-Mjj region, because of the asymmetric cuts in the transverse momenta of the two 
leading jets. 


6 Conclusions 


In this paper, we have presented the hrst complete calculation at next-to-leading order in 
perturbative QCD of one-jet and two-jet cross sections in polarized hadronic collisions. To 
this end, we extended the general formalism of ref. based upon the subtraction method 


and relevant for unpolarized incoming beams, to the case of polarized beams. The resulting 
formulae together with the matrix elements of refs. B0 
Monte Carlo code. We used the code presented in ref. 


M 


can be easily implemented in a 
for unpolarized hadron-hadron 
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collisions, and we snitably modified it in order to deal with the polarized caseQ. The program 
ontpnts the momenta of the hnal-state partons plus a weight. The momenta are eventually 
used in an analysis routine to define the physical observables. During the same computer run, 
we can obtain predictions for an arbitrary number of one-jet and two-jet quantities, using 
several jet-hnding algorithms, as well as predictions for other infrared-safe observables, such 
as shape variables. We stress that, in spite of these features, the code is not equivalent to 
a Monte Carlo parton shower generator, being the result of a fixed-order perturbative QCD 
calculation. 

Using the aforementioned code, we investigated in some detail the phenomenological 
implications of jet production at RHIC (polarized pp collisions with a maximum centre- 
of-mass energy of 500 GeV). We hrst studied the scale dependence of single-inclusive and 
double-differential jet observables, in order to determine whether the next-to-leading or¬ 
der results are reliable enough to give sensible theoretical predictions. We found that the 
next-to-leading order corrections sizeably reduce the scale dependence with respect to the 
Born results. Furthermore, the polarized next-to-leading order cross sections display a scale 
dependence comparable to that of the corresponding unpolarized cross sections. We then 
turned to the study of the dependence of our results upon the choice of the polarized parton 
densities. These quantities are very poorly constrained by the available DIS data at present, 
especially for the gluon density, which enters the dominant contributions to the jet cross 
sections in the RHIC energy range. It follows that, depending upon the specihc parton 
density set adopted, our predictions can vary by about two orders of magnitude. This fact 
in turn implies that measurements of jet observables at RHIC will be helpful in constraining 
the polarized densities. In fact, we have shown that, if the design integrated luminosity will 
be obtained, a detailed study of the polarized jet cross sections will be possible. 
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Appendix: Parton Generator for Polarized Hadron—Hadron 
Scattering 


We will briefly describe here the main differences between parton generators for unpolarized 
and polarized hadron-hadron scattering. A detailed description of the method can be found 
in ref. []^ ; here we only explicitly present the steps necessary to convert the parton generator 
of ref. 


14| into a parton generator relevant to polarized collisions. 


^The codes for both the polarized and unpolarized collisions are available upon request. 
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As already mentioned, a fake measurement function is used to disentangle the various 
infrared singularities appearing in the phase space, in order to apply the subtraction method 
to terms that contain one soft and one collinear singularity at most. In this respect, there 
is no difference between the cases of polarized and unpolarized scattering. Furthermore, 
the proof of ref. |T^, that the singularities which appear in the intermediate steps of the 


calculation eventually cancel in the sum that dehnes physically observable quantities, also 
goes unchanged. This is due to the fact that, as in the unpolarized case, the singular part 
of a matrix element has the form of a reduced matrix element convoluted with a universal 
(i.e. process-independent) kernel. 

We start from eq. 0 ; as in ref. |^, we will deal in the following with the case of iV — 1 
jet production. We write the polarized partonic cross section at NLO as 


A2) = ddSAAi; A^) + A2) + A^), 


(A.l) 


where ddiiLa is the Born contribution. All the terms in the RHS of eq. (|A-1|) are finite; 


dd^aia 2 corresponds to conhgurations with N {N — 1) partons in the hnal state. 

The iV-parton contribution reads (see eq. (A.l) of ref. H) 


7V+2 


/ 




i=3 


A 2 ) + ^ ^ 2 ) 


N+2 


\ 


(outj) I 




(A.2) 




The quantities da^^^^pfiXi, A 2 ) and dappp^j{Xi; X 2 ) can be obtained from eqs. (A.3) and (A.11) 
of ref. [0 respectively, with the formal substitution of the matrix elements squared with 


their polarized counterpart 

MW({a,}f+2) 


Ai; A2 


(A.3) 


Notice that in the unpolarized case the quantity also includes an average factor for the 

sum over the polarizations of the incoming partons, which must be removed in the present 
case. 

The (A^ — l)-parton contribution reads (see eq. (A.15) of ref. |T^) 

ddp^P~^\Xi; X 2 ) = dapp-^^\Xi; X 2 ) + dapp-^^\Xi; X 2 ). (A.4) 


The quantity ^^^(Ai; A 2 ) can be obtained from eqs. (A.16), (A.23), (A.24) and (A.25) 

of ref. 113 , with formal substitutions analogous to that of eq. (A-31) 



Ai; A 2 ), 

(A.5) 


Mpp Ai; A 2 ), 

(A.6) 


—^ M.^ns ’ Ai; A 2 ). 

(A.7) 


We stress the fact that the coefficients Q({az}]^^^) and which appear in eq. (A.16) 

of ref. [Q, are not modihed in the polarized case. The reason is the following: Q gets 
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contributions from the final-state collinear configurations, which are summed over polariza¬ 
tions, and from the soft part {z = 1) of the Altarelli-Parisi kernels relevant to initial-state 
collinear configurations, which is the same in the polarized and unpolarized cases. is 

obtained by integration of the eikonal factors, which appear in the soft limit of the iV-parton 
matrix element squared in both the polarized and unpolarized cases. The information on the 
polarizations of the incoming partons is therefore fully contained in Born matrix elements 
squared (eq. ( |A.5|) ), in the colour-linked Born matrix elements squared (eq. ( |A.6|) ), and in 
the hnite part of the virtual contribution (eq. (|A.7|) ). The colour-linked matrix elements can 
be defined exactly as in ref. using Born amplitudes summed only over the helicities of 
the final-state partons. 

The quantity ^^^(Ai;A 2 ) in the polarized case is only slightly more complicated 

than in the unpolarized case, eq. (A.26) of ref. [Q. We get 


d 

+ / <i(^O^,{0dA,a'-Z((,l-0KuK2,X2 


+ 


as 

dvr 


f rf?OfcK)<E2ff™(Jfi,Aii(l-0A-2) 


A2 Ois 
IA2I dvr 




d( AO^,(0 dA,aZ(Ki, An (1 - f 

(A.8) 


where 


^ ) *^afe(A) d(jQb(-|-, A) i d(Ta}j(^ , A), 


^2 

A 2 


cra6(A) = daabi\;+)±daab{\;-)- 


(A.9) 

(A.IO) 


The form of the operator Oab can be read from eq. (A.26) of ref. [14 




(A.ll) 


where P^{z, 0) + eP^f{z, 0) -f 0{e^) are the unpolarized Altarelli-Parisi kernels for z < 1 in 
d — 2e dimensions, and Kab define the scheme for the unpolarized parton densities (in the 
MS scheme they are equal to zero). We also have 


AOabio = eAP<,(i-e,o) 


. S6j 


Ct), 




(A.12) 
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where, as usual, 


— Pa+b+ f‘a-b+- (A. 13) 

It is now trivial to get the quantities dAaaj^a 2 which appear in eq. (|2.3|) . From eqs. (|A.1|) , 
and (POQ we get 


N+2 


( 


dAaa,a2 = + 


i=3 


+ E ‘'A'". 


N+2 


\ 


(outj) 

0.10.2,ij 


V 


J=3 

j¥=i 


+ dAfflJi'J-*”! + dAffil''!-*’-). 


aia2 


(A.14) 


The hrst four terms in the RHS of this equation can be obtained from the corresponding 
terms in the unpolarized case with the substitution 


M 


AM, 


(A.15) 


which directly follows from eqs. (|A.3|) , (|A.5| )-( |A?^) . The form for can be directly 

obtained from eq. (|A.8|). We have 




+ 


(+S 


5 ^ f <KAO,+QdAaZ((l - 

d 

Y, f <KAOaa,(QdAaf^,(K,.,(l - (,)K^), 

J 


(A.16) 


which is in fact completely analogous to eq. (A.26) of ref. |T4 


In summary, the main structure of a computer code which evaluates jet cross sections 
is identical in the unpolarized and polarized cases. No conceptual modification is required; 
in particular, there is no need to deal explicitly with the polarizations of the incoming 
hadrons/partons. The full information on the polarizations can be embedded in the ‘po¬ 
larized’ Altarelli-Parisi kernels AP, and in the ‘polarized’ matrix elements AM. Both 
quantities can be treated as black boxes, exactly like the corresponding quantities in the 
unpolarized case. 
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